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Abstract —This paper presents an analytical model for the 
fading channel correlation in general scattering environments. 
In contrast to the existing correlation models, our new approach 
treats the scattering environment as non-separable and it is 
modeled using a bi-angular power distribution. The bi-angular 
power distribution is parameterized by the mean departure 
and arrival angles, angular spreads of the univariate angular 
power distributions at the transmitter and receiver apertures, 
and a third parameter, the covariance between transmit and 
receive angles which captures the statistical interdependency 
between angular power distributions at the transmitter and 
receiver apertures. When this third parameter is zero, this new 
model reduces to the well known “Kronecker” model. Using 
the proposed model, we show that Kronecker model is a good 
approximation to the actual channel when the scattering channel 
consists of a single scattering cluster. In the presence of multiple 
remote scattering clusters we show that Kronecker model over 
estimates the performance by artificially increasing the number 
of multipaths in the channel. 

I. Introduction 

Wireless channel modelling has received much attention 
in recent years since space-time processing using multiple 
antennas is becoming one of the most promising areas for 
improvements in performance of mobile communication sys¬ 
tems [1,2], In channel modelling research, the effects of 
fading channel correlation due to insufficient antenna spacing 
and sparse scattering environments are of primary concern as 
they impact the performance of multiple-input multiple-output 
(MIMO) communication systems. 

A popular channel model that has been used in MIMO 
performance analysis is the “Kronecker” model [3-5]. In this 
model, the correlation properties of the MIMO channel are 
modeled at the transmitter and receiver separately, neglecting 
the statistical interdependency between scattering distributions 
at the transmitter and receiver antenna apertures. Measurement 
and analytical results presented in [6,7] suggest that the 
Kronecker model does not accurately model the underlying 
scattering channel, therefore it does not provide accurate 
performance results. 

In this paper, using a recently proposed spatial channel 
model [8], we develop an alternate to the Kronecker model 
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which gives channel correlation for a general class of scat¬ 
tering environments. In our proposed model, fading channel 
correlation is parameterized by the antenna configuration de¬ 
tails (spacing and the placement) both at the transmitter and 
the receiver arrays, and the joint bi-angular power distribution 
between transmitters and receivers, which models the scatter¬ 
ing environment surrounding the transmit and receive antenna 
apertures. The bi-angular power distribution is parameterized 
by the mean departure and arrival angles, angular spreads of 
the univariate angular power distributions at the transmitter 
and receiver apertures, and a third parameter, the covariance 
between transmit and receive angles which captures the sta¬ 
tistical interdependency between angular power distributions 
at the transmitter and receiver apertures. When this third 
parameter is zero, i.e., the power distribution at the transmitter 
is independent of the power distribution at the receiver, the 
proposed correlation model reduces to the Kronecker model. In 
order to model the scattering environment we propose several 
bi-angular power distributions and also find the correlation 
coefficients associated with these distributions in closed form. 
Using the proposed model, we show that Kronecker model 
is a good approximation to the actual channel when the 
scattering channel consists of a single scattering cluster. We 
also show that when the scattering channel consists of multiple 
remote scattering clusters, the Kronecker model over estimates 
the performance MIMO systems by artificially increasing the 
number of scattering clusters in the scattering channel. 

Notations: Throughout the paper, the following notations 
will be used: Bold lower (upper) letters denote vectors 
(matrices). [-] T , [•]* and [•]' denote the transpose, complex 
conjugate and conjugate transpose operations, respectively. 
The symbols S(-) and 0 denote the Dirac delta function and 
Matrix Kronecker product, respectively. The notation !?{■} 
denotes the mathematical expectation,vec(A) denotes the 
vectorization operator which stacks the columns of A, [.] 
denotes the ceiling operator and S 1 denotes the unit circle. 


II. Spatial Channel Model 

First we review the spatial channel model proposed in [8]. 
Consider a MIMO system consisting of ut transmit antennas 
located at positions x t , t = 1,2, • ■ ■ , rij relative to the 
transmitter array origin, and hr receive antennas located at 


positions y r , r = 1,2, • • • , n r relative to the receiver array 
origin. tr > max || x t || and r r > max || y r || denote 
the radius of spheres that contain all the transmitter and 
receiver antennas, respectively. We assume that scatterers are 
distributed in the far field from the transmitter and receiver 
antennas and regions containing the transmit and receive 
antennas are distinct. 

By taking into account physical aspects of scattering, the 
MIMO channel matrix H can be decomposed into determin¬ 
istic and random parts as [8] 


h = j r h s j ] t , (1) 

where Jr is the receiver configuration matrix, 

J-m r {v i) ••• Jm r {v l) 

J-m r {v 2) ••• Jm r {v 2) 

Jr = . . . , 

'J—Mr ( v ur) ■■■ ^Mr^ur) 

and Jr is the transmitter configuration matrix, 

\J — i\4r(tri) ' ' ' (^1) 

* * * \JMt (^2 ) 

Jt = . . . , 

\J— Mrp ) ’ ' ’ \JM’r i^'R'riT ) 


where J n (x ) is the spatial-to-mode function (SMF) which 
maps the antenna location to the n-th mode of the region. 
The form which the SMF takes is related to the shape of 
the scatterer-free antenna region. For a circular region in 2- 
dimensional space, the SMF is given by a Bessel function of 
the first kind [8] and for a spherical region in 3-dimensional 
space, the SMF is given by a spherical Bessel function [9]. 
For a prism-shaped region, the SMF is given by a prolate 
spheroidal function [10]. 

Here we consider the situation where the multipath is 
restricted to the azimuth plane only (2-D scattering envi¬ 
ronment), having no field components arriving at significant 
elevations. In this case, the SMF is given by 




where J n (■) is the Bessel function of integer order n, vector 
w = (|| w ||, (j) w ) in polar coordinates is the antenna location 
relative to the origin of the aperture which encloses the 
antennas, k = 2 tt/X is the wave number with A being the 
wave length and 1 = V—T. Jr is urx(2Mr + 1) and Jr 
is urx(2Mr + 1), where 2 Mr + 1 and 2 Mr + 1 are the 
number of effective 1 communication modes at the transmit and 
receive regions, respectively. Note, Mr and Mr are defined by 
the size of the regions containing all the transmit and receive 
antennas, respectively [11], In our case. 


Mr 


kerR 

A 


and Mr 


ker R 

A 


where e « 2.7183. 


1 Although there are infinite number of modes excited by an antenna array, 
there are only finite number of modes (2 M + 1) which have sufficient power 
to carry information. 


Finally, Hr is the ( 2Mr + 1 ) x (2Mr + 1) random complex 
scattering channel matrix with (£, m)-th element given by 

{H s } e m = [[ 9 (^, 

J J^x S 1 

( 2 ) 

representing the complex scattering gain between the 
(to — Mr — l)-th mode of the scatter-free transmit region 
and (l — Mr — l)-th mode of the scatter-free receiver region, 
where g(</>, <p) is the effective random complex scattering 
gain function for signals with angle-of-departure <f> from 
the scatter-free transmitter region and angle-of-arrival ip at 
the scatter-free receiver region. The reader is referred to 
[8] for more information about the channel decomposition m 


The correlation matrix of the MIMO channel H given by 
03 can be written as 

R = E {h/i t } = {J* T ® J R ) R s (Jr O Jje), 

where h = vec (H) and Rr the modal correlation matrix of 
the scattering channel, 

Rs = E j/rs/ily j , 


with hs = vec (Hr). Modal correlation matrix Rr can also 
be written as a block matrix of (2Mr + 1) x (2Mr+1) blocks, 
each of size (2 Mr + 1) x (2 Mr + 1), 


Rs = 


Rs, 1,1 
Rs, 2,1 


Rs, 1,2 
Rs, 2,2 


Rs,l, 2 M T +l 

Rs, 2,2M t +\ 


Rs ,2M r +1,1 Rs ,2Mr + 1,2 ‘ ‘ ■ Rs ,2Mr+1,2M t + 1 


where Rs,m,m' is the correlation between TO-th and m'-th 
columns of Hr. A diagonal block Rs^ m ,m gives the modal 
correlation matrix at the receiver region due to the TO-th mode 
at the transmit region whereas off diagonal blocks Rs^ m ,m' 
give the cross correlation between two distinct modal pairs at 
the transmit and receiver apertures. 


A. Modal Correlation in General Scattering Environments 

Using 0, we can define the modal correlation between 
complex scattering gains as 

jt%=E{{ H sh,m{H s };, tm ,}. (3) 

Substituting 0 in 0 gives 

e-iV-MR-Vvjtf-MR-iW&fl^i^' ( 4 ) 

where we have introduced the shorthand f 4 = ffff§i xS i- 
Assume that the scattering from one direction is independent 
of that from another direction for both the receiver and 
the transmitter apertures (WSSUS). Then the second-order 
statistics of the scattering gain function g(cf>, p) is given by 

E\g((t>,y)g*((tJ,y)} = G(^,(p)S(<j> - <f> )6(ip - ip ), 














where G{(f>,p) = E {\g(<t>, <^)| 2 } is the 2D joint azimuth 
power spectral density 2 (PSD) over departure and arrival 
angles <j> and p of the scattering channel, normalized such 
that the total scattering channel power 



pjdpdfi = 1 . 


Using this assumption, the modal correlation coefficient © 
can then be simplified to 


V’* , = rf-*- , 
im,m' lm—m' ’ 


ff G^yy-W-^e^-^dydrt), (5) 
J Js 1 xS 1 


which gives the (t, f')-th element of Rs,m,m>■ 


Since the scattering gain function <?(</>, p) is periodic in both 
(j> and p, the joint PSD function G(0. p) is also periodic 
in both <j) and p. Therefore, using the orthogonal circular 
harmonics e m(p as the basis set, p) can be expanded 
in a 2-D Fourier series as, 

^ OO OO 

G(yp) = —J 2 £ 

p= — o o q =—oo 

with coefficients 


% 



G(</>, p)e~ ipcp e^dpdfi. 


( 6 ) 


By comparing © with 0 , we can see that entries of Rs 
are given by 2-D Fourier coefficients of the joint PSD function 
G((j>, p). Furthermore, from 0 . it is evident that the entries of 
Rs are dependent on the joint PSD G{tj), p), which is usually 
parameterized by the mean departure and arrival angles <j> o, po, 
angular spreads at and a r for distributions at the transmitter 
and receiver apertures, and the covariance 


P = cov(^, p) = 


E{cj)p} - 4>opo 

(Jt (J r 


(7) 


between transmit and receive angles. 

From 0 , the correlation between Ath and f'-th modes at 
the receiver region due to the m-th mode at the transmitter 
region is given by 


= f V Rx (p)e l ( e e ' )ip dp, 
Js 1 


( 8 ) 


where V Rx {p) is the average power density of the scatterers 
surrounding the receiver region, given by the marginalized 


PSD 


Vrx{p)= f G(4>, p)d(j>. 
J s 1 


Similarly, the correlation between m-th and m'-th modes at the 
transmitter region due to the f-th mode at the receiver region 
is given by 

7m, m' = [ V Tx (ct>)e^ m - m '^dy (9) 

Js 1 


where Vtx(4>) = pjdp is the average power density 

of the scatterers surrounding the transmitter region. 


11. Kronecker Model as a Special Case 

When the covariance between departure and arrival angles 
is zero, p = 0, the joint PSD can then be expressed as 
the product of scattering distributions at the transmitter and 
receiver regions, i.e., 

G{(j>,p) =V Tx {ct>)V Rx (p). 

This separability condition leads to the well known ‘Kro¬ 
necker’ model [5,7,12], where we have 

Rs = Ft <8> F R 

with Ft the (2 Mt + 1) x (2 Mt + 1) transmitter modal 
correlation matrix and F R the (2M R +1) x (2M R + l) receiver 
modal correlation matrix. The (to, to')- th element of Ft is 
given by © and the (£, £')-th element of F R is given by ©. 

The separability of G(<j>, p) when p = 0 also yields that 

• modal correlation at the transmitter 0m,m' is independent 
of the mode selected from the receiver region, 

• modal correlation at the receiver '/' J is independent of 
the mode selected from the transmitter region and 

• correlation between two distinct modal pairs is the prod¬ 
uct of corresponding modal correlations at the transmitter 
and the receiver, 

i.i' (.y 

7 m ,m' — 7m,m'7 

In the Kronecker model, diagonal blocks Rs, m ,m and 
off-diagonal blocks Rs, m ,m' of Rs are given by F R and 
0 m ,m'F R , respectively. 

III. Bi-Angular Scattering Distributions 

In this section, we outline several examples of bi-angular 
scattering distributions along with their modal correlation 
coefficients 0 , which give the entries of the modal correlation 
matrix Rs- 


A. Uniform Limited azimuth field 

When the energy leaves uniformly to a restricted range of 
azimuth (fio — A t , </>q + A t ) at the transmitter and to arrive 
at the receiver uniformly from (po — A r ,po + A r ), then 
following Morgenstern’s family of distributions [13], we have 
the joint uniform limited azimuth scattering distribution 


Gu{<t>, v) 


i _ p{<t>-<t> o)(y- y>o) 

4 A t A r 4A?A2 


for — tf> o| < A t and \p — p$\ < A r , and 0 elsewhere. 
The parameter p is the covariance between cj> £ [— n, n) and 
p £ [—7r, 7r), which controls the flatness of Guifi, 1 ^)- In this 
case, the modal correlation coefficients 0 are given by. 


{ sine ((to - if 1 = 1' and to y to' 

sinc((£ — £')A r )e~ z ^ i '' )Vo , if and to = to' 

otherwise 


where , is given by (flOl . Note that Gu(fi,p) has 
marginal distributions Vrxifi) = 1/2A ( for <f> £ (fio — 
At,cj>o + A t ) and zero elsewhere, and V Rx (p) = l/2A r 
for p £ (po — A t ,pq + A r ) and zero elsewhere, with 


2 also called bi-angular power distribution or joint scattering distribution. 





corresponding transmit and receive modal correlation coef¬ 
ficients 7 m ^m' = sinc((m — m ')and 7^ = 
sinc((£ — £')/\ r )e~ l ^~ i '' >ipo , respectively. For A t = 7r and 
A r = 7r with p = 0 (isotropic scattering), we have uniform 
PSD Gu (</>, <p) — 1 /47 t 2 and the modal correlation coefficients 
become, 

7 m ,m' — Om-m'di-i ', 

corresponding to i.i.d. elements of FTs. In this special case, 
the correlation matrix if takes the form 

i? = (J^J?)<g> 


IV. Simulation Examples 

In this section we compare the performance of MIMO com¬ 
munication systems operating in separable (kronecker channel 
with p = 0 in <0) and non-separable scattering environments. 

We consider transmit and receive apertures of radius 0.5A, 
corresponding to 2|~7re0.5] + 1 = 11 effective communication 
modes at each aperture. Within each aperture, we place three 
antennas in a uniform circular array (UCA) configuration (3x3 
MIMO channel). The system performance is measured in 
terms of the average mutual information. Here we assume 
transmitter has no knowledge about the channel and the 
receiver has the full knowledge about the channel. In this case, 
the average mutual information is given by 


B. Truncated Gaussian Distributed Field 


A distribution that can be used to model the joint PSD is 
the truncated Gaussian bivariate distribution, defined as 


G g (<^) 


f Ig exp 


<pY 

. 2(1 — p 2 )\ ’ 


<t>,v e [-a 71- ), 


where FIq is a normalization constant such that 

/ Jgi Gg{(/>, <p)d(/>d<p = 1 and 


Q(</>, <p) 


/t -2 


2p(<t>-(l>o){<p-<po) 

(Jt CT r 


{(fi-ipo) 2 


with (j) 0 the mean AOD at the transmitter, at the standard devi¬ 
ation of the non-truncated marginalized PSD at the transmitter, 
ipo the mean AOA at the receiver, a r the standard deviation of 
the non-truncated marginalized PSD at the receiver and p is 
the covariance between receive and transmit angles, as defined 
by ®. In this case, finding modal correlation coefficients in 
closed form poses a much harder problem. However, if the 
angular spread at the both end of the channel is small, then 
a good approximation for the truncated Gaussian case can be 
obtained by integrating over the domain (—00, 00), since the 
tails of marginalized PSDs cause a very little error. Using a 
result found in [14], 


7^ m '~ ex p(*(( m - m ')0o - {£-£')vo)- 

i ( at{m — m , ) 2 — 2pata r (m — m')(£—£') + o >(£—£') 2 ))- 


C. Truncated Laplacian Distributed Field 


1 = E 



-HH 


ut 


where 7 is the average symbol energy-to-noise ratio (SNR) at 
each receiver antenna. 

Assuming FI is a positive definite matrix, a realization of 
the MIMO channel H is obtained by forming 

vec(iT) = F? 1/2 vec (W), 


where R 1 ^ 2 is the positive definite matrix square root of R 
and W is a hr x tit matrix which has zero-mean independent 
and identically distributed complex Gaussian random entries 
with unit variance. 



Similar to the truncated Gaussian distribution, an elliptical 
truncated bivariate Laplacian distribution can be defined as 

[15] 


Gl{4>,v) = lKq 


l 2Q((j),ip) 
(1 -P 2 ) 


€ [—7T, 7t), 


where FLr is a normalization constant such that 
/ Jgi Gg( 0 : (p)d<pdp = 1 and -Ko(’) is the modified 

Bessel function of the second kind of order zero. The modal 
correlation coefficients 0 for this distribution are given by 


t,e' 

^ m,m' 


exp(z(m — m)(f)Q — (£—£')ipo) 

<jt(m — m') 2 — 2 pat(7 r (m — m') (£—£') + <j%(£—£') 2 + 1 


Fig. 1. Average mutual information of 3-transmit UCA and 3-receive UCA 
MIMO system in separable (Kronecker with p = 0) and non-separable (p = 
0.8) scattering environments: bivariate truncated Gaussian azimuth field with 
mean AOD = 90°, mean AOA = 90°, transmitter angular spread crt = 10° 
and receiver angular spreads a r = {30°, 10°}. 

Figure 0 shows the average mutual information for a 
bivariate truncated Gaussian distributed azimuth field with 
p = 0.8. It was shown in [16] that performance of UCA 
antenna configuration is less sensitive to change of mean 
AOD (</>o) and mean AOA (po). Therefore, without loss of 
generality, we set (j>o = To = 90°. Also, in this simulation, we 
set transmitter angular spread a t = 10° and receiver angular 
spreads a r = {30°, 10°}. For comparison, also shown is the 
average mutual information of the 3 x 3 i.i.d. MIMO channel. 
We observe that when ay = 30°, both models give very similar 


















p 


X 


r 


m,m' 


= sinc((m - m , )A t )sinc((£ - £')A r ) + _ 

[cos((m — m')A t ) — sinc((m — m') A t )] [sinc((f - 


m — to') At A r 
- t')A r ) — cos((£ — £')A r )]. 


( 10 ) 


performance for all SNRs. When the angular spread at the 
receiver is small, e.g. ay = 10°, we can observe that the 
Kronecker model gives slightly higher performance than the 
non-separable model for higher SNRs. However, the margin 
of capacity over estimation is insignificant in comparison with 
the i.i.d. channel capacity performance. Therefore, Kronecker 
model provides a good estimation to the actual scattering 
channel when the joint scattering distribution is uni-modal. 
Reasoning for this claim will be discussed in the next section. 

A. Capacity in Multi-Modal Bivariate Azimuth Fields 

A multi-modal azimuth power distribution arises when 
there are two or more strong multipaths exist in a fading 
channel. This may be the result of multiple remote macro¬ 
scopic scattering clusters, for instance. A multi-modal bivariate 
distribution can be constructed from a mixture of uni-modal 
bivariate distributions. Fig. |3 shows a multi-modal bivariate 
Gaussian distributed azimuth field with 3 modes centered 
around (<j>o,<po) = {(-40°, 40°), (0°,-40°), (50°, 0°)}, each 
mode with angular spreads ay = ay = 5° and p = 0.8. Note 
that, in this case the effective angular spreads at the receiver 
and transmitter are larger than 5°. 


25 'l 

20 -v 


15-v 



Fig. 2. An example multi-modal bivariate Gaussian distributed azimuth field. 

We now consider the 3x3 antenna configuration setup 
discussed in the previous example. Fig. 0 shows the average 
mutual information of it applied on the multi-modal scattering 
distribution shown in Fig. |3 It is observed that Kronecker 
model tends to overestimate the average mutual information 
at high SNRs. Unlike in the uni-modal case considered pre¬ 
viously, the margin of error seen here is quite significant, 
especially at high SNRs. Following the analysis given in [7], 
we now provide reasons for Kronecker model to overestimate 
the mutual information for the scattering distribution shown 
in Fig. 0 



Fig. 3. Average mutual information of 3-transmit UCA and 3-receive UCA 
MIMO system for separable and non-separable scattering channel considered 
in Fig. [2 

The joint PSD of the Kronecker model is given by 
G(4>, <p) = V T x{f)V Rx {p), where Vtx{4>) and V Rx (tp) are 
the transmit and receive power distributions, generated by 
marginalizing G(<j>, <p). Fig. 0] shows the Kronecker model 
PSD G(4>, <p) of the scattering channel considered in Fig. Q 
Comparing Fig. ^ with Fig. |3 we can observe that G(<j>, ip) 
consist of six extra modes, corresponding to additional six 
scattering clusters. Therefore, Kronecker model introduces 
virtual scattering clusters located at the intersection of the 
actual scattering clusters. As a result, Kronecker model will 
increase the effective angular spread at the transmit and re¬ 
ceive apertures (lower modal correlation) and hence improved 
system performance. Therefore, the popular Kronecker model 
does not model the MIMO channel accurately when there exist 
multiple scattering clusters in the channel These observations 
match the measurement results published in [6]. 

Now we consider the uni-modal PSD used in our first 
simulation example. Fig.^shows the corresponding Kronecker 
Model PSD G(<t>, ip) for this channel, for ay = 10° and 
ay = 10°. In this case the Kronecker model does not introduce 
any additional virtual scattering clusters into the channel. As a 
result, no increase in the number of multipaths of the channel, 
hence both models give very similar performance. 

V. Conclusions 

We presented a MIMO channel correlation model which is 
capable of capturing antenna geometry and joint correlation 
properties of both link ends. The scattering environment 
surrounding the transmitter and receiver apertures is modeled 
using a bi-angular power distribution. We use the covariance 
between transmit and receive angles to control the joint 
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Fig. 4. Kronecker model PSD <p) = Vtx (4>YPrx (y) of the non- 
separable scattering distribution considered in Fig. |3 
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Fig. 5. Kronecker model PSD G{<j >, ip) = Vtx (<I>)'Prx (:f ) of the uni-modal 
non-separable scattering distribution used in the first example to obtain the 
results in Fig. [7] lor a r = 10°. 

correlation properties between transmit and receive angular 
power distributions. 

We showed that 2-D Fourier series coefficients of the bi- 
angular power distribution and transmit and receive antenna 
sampling points contribute to the entries of the correlation 
matrix. We proposed several bi-angular power distributions 
and their 2-D Fourier series coefficients in closed form. 

Using the proposed model, we show that Kronecker model 
is a good approximation to the actual channel when the 
scattering channel consists of a single scattering cluster. In the 
presence of multiple remote scattering clusters we show that 
Kronecker model over estimates the performance of MIMO 
systems by introducing virtual scattering clusters into the 
channel. Therefore, in this case, Kronecker model cannot be 
used to represent the channel. 
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